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precision from its vertex detector (VELO) is essential to meet the tight requirements of vertex reconstruction in B-physics. The 
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1. Introduction 

LHCb is due to commence data taking as the first 
dedicated B-physics experiment at a hadron collider 
in the second half of 2008. One of the key character- 
istics of measurements of B-meson decays is their dis- 
placed decay vertex. This requires particularly precise 
measurements of both primary and secondary vertex 
positions. For LHCb, this is performed by the Vertex 
Locator (VELO) Q]. 

With a single hit resolution of better than 10 /im [2J, 
even for tracks with minimal charge sharing, a knowl- 
edge of the VELO sensor positions on the micron level 
is vital. Studies have shown that for example a mis- 
alignment of the VELO sensors of the order of the sin- 
gle hit resolution deteriorates the precision on the B- 
meson proper time by about 50% [3]. 

A challenge to the determination of alignment con- 
stants is posed by the movement of the two halves away 
from the beam during injection. This means that more 
than once a day each VELO half is moved by 3 cm 
each way. Each time the vertex locator is moved a new 
alignment is required. 

In section 2 the techniques for alignment are ex- 
plained in detail. In the third section the performance 
of the alignment algorithms determined from a beam 
test is discussed and followed by the conclusion in sec- 
tion four. 

2. VELO alignment method 

The LHCb VELO consists of two retractable halves 
with each half having 21 modules that each support 
two semi-circular silicon strip sensors, one measuring R 
(radial) and the other $ (azimuthal) co-ordinates [2]. 
When in their nominal position, the sensors form a se- 
ries of disks, of radius 42 mm, placed normally to the 



beam axis. A cut out up to a radius of 7 mm allows 
for the beams to pass through the detector without di- 
rectly interacting with the sensors. However, the LHC 
filling procedure requires that the VELO is retracted 
by at least 30 mm from its nominal operating position. 
Thus, during injection the two detector halves are sep- 
arated by m 60 mm. 

After a precise assembly and a detailed metrology, 
the position of the active parts of the VELO during 
operation will be determined by track based software 
alignment algorithms |3J. This procedure is divided 
into three stages. First, the relative alignment of the 
two sensors on each module are determined. This is 
followed by the alignment of all modules in each of the 
two halves with respect to each other. Finally, the rel- 
ative position of the two halves is measured. 

All algorithms exploit measurements of track resid- 
uals. Here, residuals are the difference between the in- 
tercept of a fitted track with a sensor and the mea- 
sured hit position on this sensor. A systematic bias of 
the residual distribution for one sensor can then be as- 
sociated with a misalignment. 

For the relative alignment of the two sensors on each 
module it is important to note that residuals arc only 
sensitive to misalignments perpendicular to the direc- 
tion of the strip. As the strip direction varies across the 
sensor, this produces a relationship between the resid- 
ual distribution across the sensor to its displacement 
as follows: 

residuals. = —Ax cos 4>track + Ay sin 4> tr ack (1) 

residual® = Ax sin farack + Ay cos 4> track + Ajr tra $) 

with Ax and Ay the translation misalignment along 
the x and y axes respectively, and A7 the rotation mis- 
alignment around the z axis. For each sensor the resid- 
ual distribution is determined from a sample of about 
20000 tracks and then fitted according to these equa- 
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tions. The difference in the fit result gives the relative 
misalignment of the two sensors. 

Both the relative module alignment and the VELO 
half alignment use a global minimisation method based 
on the Millepede algorithm [5]. This procedure relies 
on a set of linear equations relating a track model to 
residuals, where the latter are expressed as a function 
of alignment constants (A) and their derivatives (c). 
Summing up over all tracks (tr) and dimensions of the 
track model (?) one can derive a global x 2 function: 




Etr;, 
X 3 



EtrA 



•A fc /a 2 , (3) 



with the measurements (y), the track parameters (a) 
and their derivatives (x). The minimised % 2 can then 
be written as a matrix equation. 
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V a track ) — y f(X.Y,a)J 



(4) 



Inverting the matrix V directly yields the optimal set 
of alignment constants (A). 

However, this matrix is of the size of the number of 
global parameterq_!l plus the number of tracks times 
the number of parameters per track. As one needs of 
the order of many tens of thousands tracks to align the 
modules in a VELO half, inverting a matrix of this size 
is non-trivial. The Millepede algorithm exploits the 
fact that this matrix is very sparse as its major block 
has only small blocks for each track along the diago- 
nal. Using this knowledge, the whole procedure is able 
to align the modules inside each of the VELO halves 
within a few minutes on a single CPlo- Furthermore, 
it intrinsically takes into account the correlations be- 
tween the different modules and includes the function- 
ality for constraining weak degrees of freedom, i.e. de- 
formations that are difficult or impossible to unfold 
with tracks. 

To align the two VELO halves with respect to each 
other, the same algorithm is used. However, the track 
sample used here consists of tracks that produce hits 
in both VELO halves. This is possible due to a small 
region in the x-y view (around 2% of the active area) 
where sensors from both sides of the VELO overlap. 
As these tracks produce a very stringent constraint on 
the relative position of the two halves, only a few hun- 
dred are needed to determine the alignment constants. 
These are expected to exist in sufficient quantities of 
either muon beam halo tracks or tracks from beam-gas 
interactions. As an alternative, tracks from a common 
primary vertex can be used to align the two halves. 
However, due to the additional degrees of freedom of 
the vertex fit, these only provide a weaker handle on 
the alignment parameters. 

During data taking it is foreseen to run both the 
module alignment and the VELO half alignment algo- 



1 for the modules in one VELO half: n s ; f, a ; 
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rithms after each fill when the two halves have reached 
their closed position. This process will take about 10 to 
15 minutes and will thus ensure a high precision of the 
alignment already for the first processing of the data 
which starts with a delay of a few hours with respect 
to the data taking. The relative sensor alignment is 
not expected to change during the VELO movements 
and will hence be performed at larger time intervals. 
However, all alignment constants will be continuously 
monitored. 

3. Beam test results 

A first test of the alignment methods using real data 
was performed after a beam test in late 2006. The beam 
test setup comprised ten VELO modules and used the 
same hardware and software readout chain like the fi- 
nal experiment. During different phases of this test dif- 
ferent sets of six out of the ten VELO modules were 
read out simultaneously. 

The beam used was a 180 GeV secondary pion beam 
which was directed either directly onto the sensors or 
aimed at small lead targets that produced secondary 
tracks to be detected by the VELO sensors. The targets 
were located at various positions along the position of 
the nominal beam line. 

A comprehensive overview of the performance of the 
alignment software with beam test data is given in 
Ref. |6J . A first indication of the quality of the align- 
ment algorithms can be obtained by studying resid- 
ual distributions. As described in the previous section, 
the amplitude of any observed sinusoidal distribution 
of the residuals as a function of the azimuthal angle 
(j> is directly related to the size of the translational 
misalignment of the sensors. To get a more quantita- 
tive measure, the shape of these residual distributions 
has been projected to the residual axis for all sensors, 
yielding a one-dimensional distribution whose width 
is related to the size of the remaining misalignments. 
Figure [T] shows this distribution and the width of only 
0.9 (im confirms that an alignment precision on the 
micron level has been achieved. 

In the same way, the slope of the residual distribu- 
tion as a function of the radial co-ordinate r gives the 
size of the misalignment due to a rotation around the 
z axis. Again, its projection to the residual axis can be 
related to the overall size of the remaining misalign- 
ments. Figure[2]shows the distribution of the projected 
residual distributions vs r for all $ sensors. The width 
of 2.0 /Ltm can be related via an effective radius to the 
precision on rotational misalignments, yielding a size 
of remaining misalignments of 0.1 mrad. 

Compared to the single hit resolution for the smallest 
strip pitch and for tracks with minimal charge sharing 
of roughly 9 /im these levels of alignment precision lead 
to the conclusion that no deterioration in the detector's 
performance has to be expected due to misalignments. 

One effect that can potentially limit the accuracy 
of the alignment algorithms and thus the single hit 



Quality of x/y translation | 






40 


E> I 




35 


r 




30 


E- 




25 


i. 




20 


1_ 




15 


r 




10 


E- j[ 




5 


E- 


1 S 


°2 


-, , , , I , , , , 1 , , , ,r 


tf,.r.,.UL.,t,..nt..,., 


-15 -10 


5 5 10 15 20 

projected residual in iim 



Fig. 1. Alignment quality plot showing the projection of the 
residual distribution vs <j>. The width of the distribution of 
0.9 (im shows the size of the remaining translational misalign- 
ments. 
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Fig. 2. Alignment quality plot showing the projection of the 
residual distribution vs r. The width of the distribution of 
2.0 (im shows the size of the remaining misalignments, trans- 
lating into a remaining rotational misalignment of 0.1 mrad. 

precision is asymmetric cross-talk introducing a bias 
depending on the readout order of the single strips. 
Particularly in the case of the R measuring sensors, this 
effect is strongly visible and can lead to a bias of several 
microns. It originates from the fact that the readout 
direction is such that, for some parts of the sensor, 
strips are read out towards increasing radii while, in 
other parts of the sensor, strips with smaller radius 
are read out later. A forward cross-talk introduced in 
long-distance cables, connecting the detector with the 
readout electronics, hence leads to a bias either towards 
a smaller or larger radius. 

To compensate for this, a finite impulse response fil- 
ter has been put into place. This algorithm, when per- 
formed during the digitisation of the raw data, has been 
shown to be able to correct for the cross-talk to a level 
that does no longer limit the alignment algorithms. 

A completely independent test of the quality of the 
alignment software has been performed using beam 
test data by comparing the alignment constants as ob- 
tained by the software algorithms to those measured 
by the module metrology during assembly [7] . Figured 
shows this comparison for the relative translation of 
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Fig. 3. Comparison of alignment constants determined by soft- 
ware with thus determined by a smartscope metrology. 
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Fig. 4. Reconstructed vertices as function of their position along 
the beam axis without (dashed) and with correction of mis- 
alignments. 

R and $ sensors on six modules. It shows that both 
methods are in agreement down to a level of 5 /zm 
which is comparable to the expected accuracy of the 
smartscope metrology, thus confirming the precision of 
the software alignment. 

An obvious example for the need of alignment al- 
gorithms turned out to be the reconstruction of ver- 
tices of tracks from interactions of the beam in the 
targets. Despite having two targets in a certain region 
traversed by the beam, only vertices from one target 
were reconstructed when assuming the nominal geom- 
etry of the sensors. After correcting for misalignments 
with the output of the alignment algorithms, not only 
vertices from the second target were reconstructed but 
also the vertex resolution showed a clear improvement 
(see Fig. [3J. In general, precise vertexing is vital for 
LHCb due to the tiny displacement of the B decay ver- 
tices that have to be reconstructed already at trigger 
level. 

Finally, another valuable cross-check was the study 



of the stability of the VELO alignment. Therefore, 
alignment constants have been determined while the 
detector was operated in air at atmospheric pressure 
as well as during operation in a vacuum of 10~ 3 mbar. 
As no significant change in the alignment constants 
has been observed it was confirmed that there is no siz- 
able deformation due to effects like out-gassing during 
evacuation of the detector volume. 

Another two sets of alignment constants that have 
been compared were determined from data taken be- 
fore and after the detector had been moved, similar 
to the movement during each fill of the LHC. Here as 
well, no significant change in the alignment constants 
has been seen. It may be concluded that the detector 
alignment, and hence its performance, is not compro- 
mised by the need to move the detector between fills. 

4. Conclusion 

The VELO software alignment has been presented 
with its three steps: the relative alignment of the sen- 
sors, using the shape of residual distributions, and the 
module and VELO half alignment using the global \ 2 
minimisation algorithm Millepede. With these algo- 
rithms all necessary alignment constants can be deter- 
mined within minutes on a single CPU. 

It has been shown that with a precision at the mi- 
cron level, the algorithm uncertainties are well below 
the single hit resolutions. Hence, the alignment will not 
introduce adverse systematic effects on physics observ- 
ables. An excellent agreement of software results with 
metrology data underlines the precision in the deter- 
mination of the alignment constants. Finally, a good 
stability of the VELO modules under movements and 
evacuation of the detector volume has been demon- 
strated. 
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